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Abstract 

Extending the Kruppa's prescription for the continuum level density, we have recently im- 
proved the BCS method with seniority- type pairing force in such a way that the effects of 
discretized unbound states are properly taken into account for finite depth single-particle po- 
tentials. In this paper, it is further shown, by employing the Woods-Saxon potential, that the 
calculation of spatial observables like nuclear radius converges as increasing the basis size in the 
harmonic oscillator expansion. Namely the disastrous problem of a "particle gas" surrounding 
nucleus in the BCS treatment can be circumvented. In spite of its simplicity, the new treatment 
gives similar results to those by more elaborate Hartree-Fock-Bogoliubov calculations; e.g., it 
even mimics the pairing anti-halo effect. The obtained results as well as the reason of con- 
vergence in the new method are investigated by a variant of the Thomas- Fermi approximation 
within the limited phase space which corresponds to the harmonic oscillator basis truncation. 
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I. INTRODUCTION 



The advent of new radioactive beam facilities makes it increasingly interesting to study 
unstable nuclei near the neutron or proton drip line. One of their specific features is weak 
binding of constituent nucleons, which leads to spatially extended nuclear profiles and a 
striking phenomenon of the neutron halo Among many issues expected in researches 



of nuclei far from the stability js, 3|, the pairing correlation plays a special role because 
virtual scatterings into continuum states occur more easily. The basic quantity concerning 
the spatial distribution of nucleons is the nuclear radius, which is also a prerequisite for the 
analysis of reaction cross sections However, its theoretical evaluation in the presence 
of pairing correlation is not straightforward due to the fact that the continuum states, 
into which weakly bound nucleons virtually scatter, occupy infinite volume. In fact, the 
calculated radius is unreasonably large or divergent because of finite occupation probabil- 
ities of unbound states in the simple-minded BCS treatment; the so-called "neutron gas" 
problem, the solution of which requires the more sophisticated Hartree-Fock-Bogoliubov 
(HFB) theory [iiel. 

A similar problem exists for the Strutinsky shell correction calculation with finite depth 
potential [7] , where the smoothed part of binding energy does not converge when increas- 
ing the size of the single-particle basis: The finite occupation probabilities of continuum 
states are required for extracting the smooth part, but the level density of unbound states 
is infinite, which leads to divergence of the Strutinsky smoothed energy. An efficient 
method to avoid this problem was introduced by Kruppa 8|], and was used in the shell 
correction method 9!]. The idea of the method is to calculate the so-called continuum 
level density as a difference between the level densities obtained by diagonalizing the full 
Hamiltonian including a finite depth potential and the free Hamiltonian. By employing 
the harmonic oscillator expansion, it was shown for the Woods-Saxon potential that the 
smoothed energy calculated by the Kruppa method well converges as increasing the size 



of basis {9]. 

Recently, three of the authors of this paper have proposed a new method of BCS 
calculation [10], which is free from the divergence as increasing the model space, by 



utilizing the similar idea to Kruppa's. In Ref. 



10| not only the treatment of pairing 



correlation but all the other procedures in the microscopic-macroscopic method for the 



calculation of nuclear binding energy are reexamined and improved. In this paper, we 
further show that the new method is capable of calculating the spatial observables like 
the nuclear radius without the problem of particle gas surrounding nucleus. 



II. KRUPPA PRESCRIPTION FOR EXPECTATION VALUES OF ONE-BODY 
OPERATORS 



Basic idea 



In the Kruppa's prescription 



a 



the level density is replaced as 

M M 



9{e) 9''{e)=g{e)-goie)^Y.^{e-e,)-J2S{e-e'^), (1) 

i=i j=i 



where and are the eigenvalues of the full and the free Hamiltonians, respectively. 
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In the following, we are mainly concerned with neutrons, but it should be reminded for 
protons that the Coulomb potential is included in the free Hamiltonian, i.e., "free" here 
means that the nuclear potential is left out. The eigenvalues in Eq. ([T]) are calculated by 
diagonalizations with the harmonic oscillator basis, and M is the total number of the basis 
states commonly used in the two diagonalizations. Both the full and free level densities, 
g{e) and go{€), are divergent as increasing the basis size (M — > oo) for the positive single- 
particle energy, e > 0. It is shown for the Woods-Saxon potential [8] that their difference 
g^{e) remains finite, and converges, in the spherical case, to the well-known expression in 
terms of the scattering phase shift 6ij{e), 

71 ^-^ de 

ij 

The great merit of Eq. is that it can be easily applied to the deformed cases (see 
Ref. [8| for the corresponding expression to Eq. ([2]) in terms of the scattering S-matrix in 
such general cases). 

Extending the idea of subtracting the free contribution in Eq. ([T]), we propose to calcu- 
late the expectation value of an arbitrary one-body observable O by a similar replacement 

as 

(O) (0)K = (O) - (0)o, (3) 
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where the first (second) term is a summation of the expectation values with respect to 
the wave functions calculated by diagonalizing the full (free) Hamiltonian multiplied with 
the occupation numbers. Note that, in the independent particle approximation, e.g., the 
Hartree-Fock (HF) theory, the second term does not contribute for bound systems, i.e., 
if the Fermi energy is below the particle threshold. The free contribution manifests itself 
when the occupation probabilities of continuum states are non-zero due to the residual 
interactions, e.g., in the case of BCS treatment for pairing correlations. 

Although generalizations are possible, we consider in this work the following seniority- 
type separable pairing interaction. 



pair 



(4) 



i>0 



where i represents the time- reversed state of i, forming a time reversal pair (ii) (they are 
degenerate in energy, ej = e^), and X]j>o i^^ans the sum is taken over all the pairs. A 
cutoff of the model space is necessary for such a pairing interaction, and it is realized 
by introducing a smooth cutoff function ll|], /c(e) (we use a different form from that of 
Ref. [Ill), defined by 



1 + erf 



e - A + Ai 



d. 



cut 



1/2 



erf 



-e + A + Au 



d, 



cut 



1/2 



(5) 



where erf(x) 



e * (it is the error function. We use the cutoff parameters of the 



pairing model space, Au = Ai = 1.2 hu and dcut = 0.2 hu, with hui = 41/ A^^^ MeV. The 
predefined parameter A in the cutoff function ([5]) is chosen to be the smoothed Fermi 
energy obtained by the Strutinsky smoothing procedure, see Ref. 10|] for details. 

The gap equation can be derived from the variational principle with the number con- 
straint, 

A = G{P^), and N = {N), (6) 

where is the neutron or proton number to be fixed. Applying the prescription ([3]), the 
modified gap and number constraint equations are obtained: 



2 
G 

N 



E 



j>0 



/c(e.)^ /c(e°)2 



(7) 
(8) 



i>0 



Here the BCS quasiparticle energy and occupation probability are given, as usual, by 

E{e) = v/(6-A)2 + /,(e)W, v'{e) = ^1 - ^) • (9) 

From these equations, the pairing gap and the chemical potential, (A, A), are determined 
for a given pairing strength G. The second terms in the square brackets in Eqs. 
and ([8]) are the extra (negative) contributions form the free spectra. We call this new 
way of the BCS treatment Kruppa-BCS method, and the resultant equation Kruppa-BCS 
equation. The consequences of this Kruppa-BCS equation have been investigated in detail 
in Ref. jl^. As long as the chemical potential is well below the free particle threshold 
(minimum of {e^,i = 1,2, ...}), the Kruppa-BCS equation has a unique solution and can 
be solved in the same way as the ordinary BCS equation. It is shown that the calculated 
pairing gaps with this new method converge to reasonable values for large basis size in 
contrast to those with the ordinary BCS equation, which strongly depend on the ^size of 
basis and hardly applicable to the calculations with large basis size 
a comprehensive discussion. 



12|. See Ref. [10( for 



B. Root mean square radii and deformation parameters 

Once the pairing gap and the chemical potential, (A, A), are obtained by the Kruppa- 
BCS equation ([7]) and ([8]), the expectation value of one-body observables can be calculated 
within the BCS treatment; for example, the root mean square (rms) radius by 

{rY = {r') - {r\ = ^ ^ [(^|r^K) 2v'{e,) - {z\r'\t)o 2v\e'^)] , (10) 

where the diagonal matrix elements (i|r^|i) and {i\r'^\i)o are with respect to the Woods- 
Saxon and free spectra, respectively. Figure [T] shows the calculated neutron radii for 
a stable nucleus ^^°Sn and an unstable nucleus ^^''Sn in the lower and upper panels, 
respectively. They are depicted as functions of the maximum value of the harmonic 
oscillator quantum number, iV™^"", specifying the model space. As for the Woods-Saxon 
potential, we use the parameter set recently developed by Ramon Wyss [13], which gives 
similar density distributions, both for neutrons and protons, to those obtained by Skyrme 
and Gogny HF calculations. In order to see the impact of the pairing correlation on the 
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FIG. 1: Neutron's rms radii as functions of the maximum harmonic oscihator quantum number 
N^^^, which specifies the size of basis. Those calculated by the ordinary prescription, (r^)n, and 
the Kruppa prescription, (r^)^ in Eq. (|10p . are included for a near drip-line nucleus ^^'^Sn, the 
panel (a), and for a /3-stable nucleus ^^'^Sn, (b). Both nuclei are spherical. The neutron pairing 
gap is fixed to be either = 0.5, 1.0, or 1.5 MeV. 



radius, the pairing gap in this calculation is kept constant to the values An = 0.5, 1.0, 
and 1.5 MeV, and only the number constraint equation ([8]) is solved for determining the 
chemical potential An. The Kruppa mean value (r"^)^ in Eq. f llOp and the ordinary mean 
value (r^)n are compared. Note that for the calculation of the latter, the ordinary number 
equation is solved so that the chemical potentials in the calculations of two radii, (r^)^ 
and (r^)n, are generally different. As is quite evident from Fig. [H the calculated radii by 
the ordinary method diverge as increasing the basis size N^^^. The divergence is more 
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rapid for larger pairing gaps. There is no way to obtain meaningful results even for the 
stable nucleus ^^^Sn. This is nothing but the problem of the "neutron gas" surrounding 



the nucleus 



Ml- 



In contrast, the radii calculated by the Kruppa method converge to 



definite values. Comparing with Fig. 15 in Ref . Q , our Kruppa results nicely corresponds 
to those of the HFB calculation, although the coordinate space is utilized in Ref. and 
the box size instead of A^™^^ is changed to demonstrate the convergence. 

PR 

It is sometimes recommended to be content with a small model space [3, 115|, like 
iV™^ ~ 12, for finite depth potentials, because the plateau condition for the shell cor- 
rection energy is not met with larger model spaces. Such a backward idea may work for 
stable nuclei, but it is clear from Fig. [1] (a) that the space N^^^ ^ 12 is not enough 
to calculate accurately the rms radius of ^^°Sn, and using a larger model space requires 
inevitably the Kruppa method. 

Another interesting feature of the Kruppa calculations seen in Fig. [1] is that the radius 
of stable nucleus stretches for larger pairing gaps while that of the nucleus near drip-line 
shrinks, although the absolute amount of changes is very small. The stretching of the 
radius due to the pairing correlation is well-known: The pairing induces the couplings to 
particle states above the Fermi surface, whose rms radii are larger in average than those of 
hole states, and increased occupation probabilities of particles and decreased occupation 
probabilities of holes make the expectation value of radius larger (see Fig. H] (c) and 
discussions in Sec JIID|) . In spite of this effect, the radius of near drip-line nucleus ^^°Sn 
shrinks. We believe that this is related to the interesting "pairing anti-halo" effect [l^ in 
the HFB theory, and will come back to this point later in Sec. Ill E[ 

In order to further illustrate the importance of the Kruppa prescription of Eq. (jS]), 
presented in Fig. [2] are examples of calculation for the quadrupole deformation parameter 
defined as a ratio of the quadrupole moment to the mean square radius, 

_ An {r^Y,o) _ An (r^F^o)^ _ 47r (r^F^o) - {r'Y,o)o 

In the upper and lower panels ((a) and (b)) the results for ^^^Dy and '^^'^Dy are depicted, 
respectively, as functions of the basis size. Both nuclei are axially deformed, and the de- 
formation parameters used in the Woods-Saxon potential are {P2, Pa) = (0.244, 0.053) and 
(0.269,0.035), respectively, which are obtained by the improved microscopic-macroscopic 
method of Ref. [lOj. As in Fig. [H the fixed pairing gaps of An = 0.5, 1.0, and 1.5 MeV 
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FIG. 2: Neutron's deformation parameter (jlip as functions of the maximum harmonic oscihator 
quantum number N^^^. Those calculated by the ordinary and the Kruppa prescriptions are 
included for a near drip-line nucleus ^"^^Dy, the panel (a), and a /3-stable nucleus ^^"^Dy, (b). 
Both nuclei are axially deformed. The pairing gap, = 0.5, 1.0, and 1.5 MeV, are employed. 

are used. Apparently, the quadrupole moment as well as the radius diverges as increas- 
ing the basis size due to the neutron gas surrounding the nucleus, if being calculated by 
the ordinary BCS treatment. One might expect that their ratio would converge to the 
correct value. It is not the case, however, even in stable nuclei, and the error is larger 
in neutron rich nuclei, amounting to about 10% for An = 1.5 MeV and N^^^ = 30 in 
^°^Dy. If the Kruppa prescription is used, the results with N^^^ ^ 12 are already suf- 
ficiently accurate. Namely, compared with the radius or the quadrupole moment itself, 
the convergence of the deformation parameter as their ratio against the increase of the 
model space is much faster. Therefore it is better to use the Kruppa method even for 
calculations with A^™^"" ^ 12. 
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C. Semiclassical consideration 



In the case of the single-particle level density, the convergence of the Kruppa density ([Tj) 
to the exact density ([2]) in the limit of infinite model space, M — > oo (continuum limit), can 
be proved rigorously a, [l3,[l8|. Recently, we have shown the convergence more pictorially 
by employing a variant of the semiclassical (Thomas- Fermi) approximation. We call it the 



oscillator-basis Thomas- Fermi (OBTF) approximation 10|], which is suitable to treat the 
problem with the truncation in terms of the harmonic oscillator basis. In the following, we 
briefly sketch this OBTF approximation, and show that the Kruppa prescription ([3]) for 
the expectation value of any spatial observables like the radius is convergent as increasing 
the basis size, A^^c"" ^ ^• 

The basic quantity in the Thomas-Fermi approximation is the phase space distribution 
function [191 . Because the momentum distribution is isotropic in the nuclear ground 
state [19|, the distribution function depends only on the magnitude of the momentum 



variable, which one can transform to the energy e. Thus the distribution function for the 
single-pg 

0{x), as 



p2 

single-particle Hamiltonian, H = h V(r), is written, using the Heaviside step function 

2m 



UAr, e) = ^-^J^ \e - V{r)\'/' 6 (e - V{r)) , (12) 



with which the Thomas- Fermi level density can be obtained as g^^ (e) = J f^^ (r, e)d r. In 
the OBTF approximation, the distribution function is modified to incorporate the effect 
of the limited phase space corresponding to the truncated harmonic oscillator (HO), i.e., 

/oB (r, e) = (r, e) 6 (e^fl^) + V{r)-e). (13) 

Here the quantity ej^^^(r) is the local maximum kinetic energy of the isotropic HO poten- 
tial, Vho('") = ^mwV^ 

= l^^' {RL. - r') e (i?^,. - r) , (14) 

where the maximum radius -Rmax allowed in the HO potential is specified by the cutoff 
energy ecut, related to the maximum HO quantum number N'^f^ of the basis truncation; 



i^max = l/^, e,,t = huj [{NZT + l)(iVo^r + 2)(iV-r + 3)]'/' . (15) 



(Although the anisotropic HO basis is utihzed in the actual calculation, it is enough 
to consider the isotropic HO potential for proving the convergence.) Thus, the infinite 
model space limit N^^^ — > oo is realized by -Rmax — > c>o in fosir, e), and the OBTF level 
density g^^^ie) = J /ob(^5 e)(fr is finite and calculable as long as -Rmax (or N^^^) is kept 
finite. Note that this is not the case for the usual Thomas- Fermi quantities; g^p{e) and 
/pp(r, e) are both infinite for the energy e above the particle threshold due to the infinite 
spatial volume of the phase space. In the same way, the OBTF distribution function for 
the free Hamiltonian, f^^{r,e), is defined by dropping the potential V{r) (or replacing 
it with the Coulomb potential for protons). The Kruppa level density in the OBTF 
approximation, S'obI^) = 9ob{^) — A'obI^)' is obtained by the spatial integral of f^j^ir, e) = 
foB^r, e)— /oB(r, e). Both level densities, (^qbI^) S'ob^^)' diverges as -Rmax — ^ oo, but the 
Kruppa density fi'^g(e) remains finite for the potential V{r) that vanishes rapidly enough 
as |r| — i- oo like the Woods-Saxon potential. It has been shown lOj that, for single-particle 



energies satisfying < e < e^^l^{r) + V{r) (assuming V{r) < 0) everywhere, 



{e-V{r)f^^ -e^/^] d^r, (16) 



asymptotically in the limit of -Rmax oo. 

It is straightforward to apply the OBTF approximation to the calculation of physical 



observables with pairing correlation (see, e.g., Ref. [19|] for the semiclassical treatment of 
the BCS theory). For a spatial one-body observable, 0{r), its Kruppa OBTF expectation 
value can be expressed as 



oo 



(0)^B= / 9^,m;e)2v\e)de, (17) 

J — oo 

where the BCS occupation probability f^(e) is given in Eq. ([9]), and the distribution of 
(O) with respect to the energy is defined by 

9oM0);^)^ J 0{r)C{r,e)d\ (18) 

which reduces to the level density fi'Qg(e) for the case = 1. The problem of "particle 
gas" occurs at energies above the particle threshold, e > 0, as in the same way as the 
level density. Then, similarly to Eq. (fT6!) . the quantity (fT8|) converges asymptotically as 

-Rmax ^ oo to 

\3/2 



So\((0);^)«^/[(.-vW)"^-."= 



0{r)d'r, (19) 
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for the spatial observable 0{r) like the radius or the quadrupole moment (assuming the 
upper cutoff energy of the pairing model space by the function ([5]), ~ A+A„, is smaller than 
e[^^^(r) + V{r) everywhere). This guarantees the convergence of the Kruppa expectation 
value (0)03 in the BCS treatment. 



D. Convergence in the Kruppa prescription 



In order to see how the calculated rms radii converge as functions of N^^^, we have 



done many test calculations [20|] . Generally, the rate of convergence is fast for stable nuclei 
as is already shown in Fig.[T](b), while it is slow for neutron rich nuclei depending on what 
kind of orbits exist near the Fermi surface. Figure [3] includes four selected examples of 
nuclei near the neutron drip-line (note the extension of the abscissa to N^^^ = 40 and the 
enlarged scale of ordinate compared with Fig. [1]). All nuclei except ^°^Dy are spherical. 
The radius has converged within accuracy of about 0.5% already at N^^^ ^ 20 in most 
cases where the pairing correlation is effective (A > 1 MeV). The very slow convergence 
in the weak pairing correlations, e.g., the case of An = 0.5 MeV in ^^Ni (panel (a)) is 
due to the filling of the "halo" orbit near the Fermi surface, e.g., z/3si/2 (and v2d^i2) in 
^^Ni. The description of such spatially extended wave functions requires many oscillator 
basis states. Another striking feature shown in Fig. |3] is that the calculated radii with 
the largest basis size [N^^^ = 40) are always smaller for larger pairing gaps: It is found 
that this calculated feature is rather general for nuclei near the drip-line, and mimics 
the pairing anti-halo effect in the HFB theory. We discuss this point in more detail in 
Sec. [mi 

In Ref. 10| the OBTF level density is compared with the Strutinsky smoothed level 
density to clarify the convergence mechanism by the Kruppa prescription, i.e., by sub- 
tracting the contributions of free spectra. It is instructive to perform a similar analysis 
for the mean square radius. The microscopic quantity corresponding to the distribution 
of (O) = (r^) in Eq. (JUD is 

9''{{r'y,e) = J2 W\r'\^) ' e.) - {t\r'\t)o 5(e - 6°)] , (20) 

i>0 
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FIG. 3: Convergence of neutron's rms radii calculated by the Kruppa method as functions of 
A^™^"" for various nuclei near the drip-hne, ^^Ni (a), ^^^Zr (b), 208Dy (c), and ^ssjjg (^^^y Only 
208]3y axially deformed. The values of neutron pairing gap used are An = 0.5, 1.0, 1.5, and 
2.0 MeV. 



with which the mean square radius in Eq. f|TO|) is calculated by integral 



(21) 



In Fig. m we depict the Strutinsky smoothed quantity of Eq. (!20l) and its OBTF ap- 
proximation (panel (a)), as well as the corresponding level densities (panel (b)). The 
large model space with N^^^ = 40 is used. As for the OBTF approximation, we 
use Eq. (fT9|) with O = for e > and the conventional Thomas-Fermi expres- 



V{r)) 



1/2 



e^/^ in the square brackets in Eq. f|T9|) being replaced with 



sion. I.e., (e 

|e — V{r)\^^'^ 9(e — V{r)), for e < 0, considerin g^ on ly the central part of the Woods-Saxon 
potential with spherical symmetry for V{r) [10]. The sharp cusp behaviors at e = 
seen in all panels of Fig. H] are due to the threshold effect characteristic in the OBTF 
approximation 



As it was already discussed in detail in Ref. 



10|, apart from the oscillation and the 
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FIG. 4: The distribution of (r^) with respect to the single-particle energy e (c.f. Eqs. (fTSj) . (f20]) ). 
panel (a), the level density, panel (b), and the rms radius of orbits at e defined in Eq. (p2|) . 
panel (c), as functions of e for the spherical drip-line nucleus ^^®Zr. The OBTF and Strutinsky 
smoothed ones with the smoothing parameter 7 = 0.8^ and 1.8 huj and with the six order of 



the curvature correction polynomial (p = 3 in Ref. 
basis size. 



are included. A'^^'' 



40 is used for the 



cusp at e = 0, the OBTF approximation reproduces nicely the microscopic level density 
smoothed with 7 = 1.8^. It also reproduces, but less nicely, the quantity (?((r^);e). 
The result of the smoothing with 7 = 1.8/ia; is larger than the OBTF approximation in 
—20 MeV < e < — 2 MeV (where the error becomes largest) and e > 2 MeV (from where 
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there is little contribution because of the vanishing occupation probability f^(e)). It is 
helpful to introduce a physically more meaningful quantity, "rms radius of orbits at the 
single-particle energy e", defined by 

(r2)V2(e) ^ v/p^(pR7?^, (22) 

which is also depicted in Fig. H] (c). The corresponding quantity in the HFB theory is the 
rms radius of the canonical basis states plotted versus the expectation value of the mean 
field Hamiltonian, whose gradually increasing tendency above the particle threshold is 



shown in Ref. 2l|. From Eqs. flTBl) and flTI?]) (and corresponding ones in the conventional 
Thomas- Fermi approximation for e < 0), the Kruppa OBTF radius (r^)^/^(e) in Eq. fl22l) 
monotonically increases in e < 0, takes a maximum value ('"^)^ at e = and goes to an 
asymptotic value (r^)y^ as e — t- oo, where (assuming the spherical symmetry) 



(r^)^ = j ^-V{r)r^dr/ j ^-V{r) r^dr, (23) 
{r^)v = I {-V{r))r^dr/ j {-V{r))r'^dr. (24) 



The microscopic Strutinsky smoothed radii in Fig. H] (c) behave similarly to the OBTF 
expression except a complete smearing out of the cusp (with 7 = l.'ihuj) and an increasing 
overestimation at positive energies. The monotonically increasing trend of the quantity 
(r^)^/^(e) in e < (smoothing width 7 = 0.8 /tw) clearly indicates that the orbit with 
larger energy has larger radius in average, which leads to the increase of the total nuclear 
radius for larger pairing gap. 

The analyses in this subsection have shown how finite results for rms radii are obtained 
with the Kruppa prescription in the OBTF approximation and in which ways the quantum 
mechanical results deviate from the results of this approximation. 



E. Pairing anti-halo like effect 

It is well-known that the occupation of the weakly bound orbit with low orbital an- 
gular momentum, typically s-orbits, leads to the nuclear halo and many examples are 
experimentally observed in light nuclei However, it is speculated that the increased 
pairing correlation binds these halo orbits more tightly and prevents the appearance of 
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halo phenomena in heavier nuclei [16[|. To confirm this speculation is one of the most 
interesting subjects both theoretically and experimentally. 
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FIG. 5: Neutron's rms radii calculated by the Kruppa method with N^^^ = 40 for the Ni- 
isotopes near the drip-line. The Kruppa-BCS equation is solved to obtain the selfconsistent 
jap, where the pairing force strength is determined by the improved average gap method 



pairing 
of Ref. 



10|. The results without including the pairing correlation are also included. 



In this respect, the Ni-isotopes are very interesting, whose drip- line may possibly be 
reached in near future. The coordinate-space HFB calculation in Ref. 221, as an example, 
showed that indeed an abrupt increase of radii in the HP calculation is diminished by 
including the pairing correlation. 

We have calculated the same quantity with the Kruppa method in the following way. 
The calculated neutron drip-line based on the original Woods-Saxon potential 13| is 
slightly inside of the expected position. Therefore, the modification of its depth has 
been done according to a new systematic method developed in Ref. lOj; then the drip- 
line isotope is the A = 92 nucleus (the A = 94 nucleus is unbound), while it is the 



A = 90 nucleus in the calculation of Ref. 



22| . The Kruppa-BCS equation is solved for 



each nucleus, and the full microscopic-macroscopic (Strutinsky shell correction) method 
is applied with the results that all the isotopes are spherical. The average pairing gap 
method with A = 13/\/A MeV is used to fix the strength G, and the calculated pairing 
gaps are ^ 1.2 to 1.7 MeV, monotonically increasing from the A = 80 to A = 92 
isotopes. We depict the result of our Kruppa calculations for radii of Ni-isotopes in 
Fig. |5l Although the absolute values of radii are slightly larger in our calculations, the 
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behavior of how the radii increase is very similar to the calculation shown in Fig. 7 of 



Ref. 



22|. 



It may be worthwhile mentioning that a BCS treatment with including only the reso- 
nance orbits (the resonance-BCS) is tested in Ref. [22i]: The calculated radii are not good 
approximation to those of HFB (no pairing anti-halo effect is reproduced), although the 
obtained binding energy is very good. This clearly indicates the subtlety of the neutron 
gas problem in the conventional BCS treatment, and we should be very careful to calculate 
such quantities for weakly bound systems 



The pairing anti-halo effect in the HFB theory [16[ is a result of real shrinkage of the 
hole component of the quasiparticle orbits due to the selfconsistent modification of the 
potentials. Apparently, such effect is absent in the Kruppa prescription because the BCS 
treatment of quasiparticle does not change the spatial distribution of each single-particle 
orbit. Then it is surprising that a similar shrinkage of the rms radius comes out from the 
Kruppa-BCS calculation. The key to understand the reason is the rms radius (r^)^/^(e) 
defined in Eq. (!22l) and depicted in Fig. H] (c). Note that the total mean square radius is 
calculated as an integral ( l2Ti) with the weight of the BCS occupation probability f^(e), 
which is a step function at e = A smeared with the pairing gap A. Therefore, if the 
quantity (r^)^/^(e) is larger for e < A and smaller for e > A, the increase of A leads to the 
effective shrinkage of the total mean square radius. In fact, this behavior, i.e., the local 
decrease near the particle threshold e ~ 0, is observed in the result of smaller smoothing 
width = 0.8 huj in Fig. g] (c). 

We show in Fig. [6] the same quantities for ^^Ni and ^^'^Dy, but are enlarged to see 
them more closely in the region near the particle threshold, which is most influential for 
pairing correlation in drip-line nuclei. To show further details, the Strutinsky smoothed 
results with the smaller averaging parameter 7 = 0.3 hw ^ 2A are also included. As 
is depicted in Fig. |4] (c) and Fig. [6l we have found it quite general that the rms radius 

is a decreasing function near the particle threshold e ^ 0. This 
trend comes from two factors: One is that, just below the threshold, the rms radius is 
larger than the classical OBTF values due to the quantum mechanical effect of the weak 
binding. Another is that the rms radii of the free spectra are larger than those of full 
spectra just above the threshold, and the subtraction leads a considerably smaller total 
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FIG. 6: The same as Fig. [J] (c) but in an enlarged scale, and for other unstable nuclei, '^^Ni (a) 
and ^'^^Dy (b). The Strutinsky smoothed results with the smoothing parameter 7 = 0.3^ and 
0.8 hoj are included. 

radius than the OBTF values given in Eq. ( 12^ . The former factor is natural but we do 
not understand the physical reason of the latter factor. Note that even if the (bound) 
halo orbit exists above the Fermi level, A < e(halo) < 0, the subtraction effect of the free 
spectra is larger, and consequently the total radius shrinks near the drip line, —A < A. 
Thus, the subtraction of the free contributions is essential for the reduction of radius in 
the Kruppa method. 

The actual amount of decrease strongly depends on the nature of orbits just below the 
threshold; if the orbits are of halo nature, the amount is larger so that the pairing shrinks 
the rms radius more strongly. As is shown in Fig. [3l the rms radius of ^^Ni having a halo 
orbit z^3si/2 more strongly shrinks than that of ^°'^Dy when the pairing is increased, which 
clearly corresponds to the stronger decrease of the quantity (r^)^/^(e) shown in Fig. |6l 
the deformation in ^°^Dy also prevents the appearance of strong halo orbits due to the 
mixing of orbital angular momenta. It is now clear that the mechanism of the shrinking 
rms radius when increasing the pairing correlation is quite different in our Kruppa-BCS 
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method from that in the HFB theory. From the discussion above, however, the essential 
ingredient for the conspicuous shrinking effect is the presence of the halo-hke extended 
orbits below the Fermi energy, which is common to the situation where the strong pairing 
anti-halo effect is expected in the HFB approach. In this way, the Kruppa-BCS calculation 
well mimics the pairing anti-halo effect. 



F. Neutron skin 



Since the radius can be accurately calculated by employing the Kruppa prescription, it 
seems also meaningful to apply it to the calculation of the neutron skin. We compare the 
calculated neutron skin thickness of the Sn-isotopes with recent experimental data 23 1 
in Fig. [71 Here the neutron skin thickness Arnp is a difference of the neutron and proton 
rms radii, 

Arnp = {tYJ' - {r^)y\ (25) 

and they are calculated by the Kruppa method with N^f^ = 20, which is large enough 
for the (not very neutron rich) isotopes shown in the figure. 
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FIG. 7: Neutron skin thickness calculated by the Kruppa method with N^^^ = 20 for the 
Sn-isotopes. Experimental data {23^] are also included. The same calculational procedure is 
used as that in Fig. [5l 



The procedure of the calculation is the same as in Fig. [5] for the Ni-isotopes, see 
Ref. [10| for details. It is interesting that the experimental data are nicely reproduced 
although the Woods-Saxon potential used in this work [l^ is not particularly aimed to fit 
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the quantity like the neutron skin thickness. (However, there are considerable ambiguities 
in the experimental neutron skin thickness depending on which types of experiments are 
used to extract it 124, 



25|.) For the Sn-isotopes, a recent Skyrme HF+BCS calculation 



with the SLy4 functional gives slightly smaller neutron skin thickness, see, e.g., Ref. 26 |. 
It is pointed out that the neutron skin has strong correlations with the coefficient of the 



asymmetry energy in the mass 
nuclear matter, see, e.g., Refs. 



'ormula and with the equation of state of the asymmetric 
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The determination of skin thickness is therefore 
very important both theoretically and experimentally. We hope that one can utilize the 
Kruppa method used in this work as an equal alternative for such an investigation if one 
optimizes the parameters of the Woods-Saxon potential as well as the liquid drop model. 

G. Other observables 



As an example of different kinds of observables, we consider the moment of inertia about 
the X-axis (perpendicular to the symmetry axis) employing the Kruppa prescription. 



J'^ = J- Jo- 



(26) 



Here the free contribution, for example, is given as 



Jo 



i,j>0 



E{4) + E{e^) 



(27) 



with {i\Jx\i)o being the matrix element of the operator with respect to the free spectra, 
and u{e) = a/1 — f (e)^. Eqs. (126!) and ( p7|) are obtained by the cranking method with 
taking the limit of zero rotational frequency: 



J 



K 



lim 



lim 



(2^ 



l^rot '^rot 

where the state Iwrot) is the cranked mean field associated with H' = H — ujmtJx, and 
||^rot)o is the same but for the free Hamiltonian. Thus, the moment of inertia is essentially 
the expectation value of the angular momentum operator, which is composed not only of 
the coordinate but also of the momentum variables. 

The effect of weak binding of constituent nucleons on the moment of inertia in neutron 



rich nuclei has been investigated in Ref. 



291 ] by the coordinate-space HFB method (the 
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FIG. 8: Neutron's contributions to moment of inertia as functions of the neutron pairing gap 
calculated by the Kruppa and ordinary methods for the drip-line nucleus ^"Mg. The size of 
model space is specified by = 30, with which the results well converge. In the same way 



as Ref. [29[, are included two types of calculations, one for a weakly bound situation. An ~ —0.3 
MeV, and another for a deeply bound situation. An ~ —9.0 MeV. The original depth of the 
Woods-Saxon potential is modified in each case. See text for details. 

particle-hole channel is approximated by using the Woods-Saxon potential). It has been 
reported that the dependence of moment of inertia on the pairing gap is much stronger in 
drip-line nuclei, e.g., in ^°Mg. We have done similar calculations but by using the simple 
BCS with the Kruppa prescription f l26p and without it. The results of neutron moment 
of inertia, J'n, are depicted in Fig. [HI The used deformation parameters are = 0.3 
(/34 = 0.0) 29|. Two calculations are presented in Ref. I29I] with the neutron chemical 
potential being An ~ —0.3 and —9.0 MeV at An ~ 0, the latter of which corresponds 
to a stable nucleus and is artificially realized by deepening the potential. We have done 
the same adjustment to match with these calculations; the depth of the original Woods- 
Saxon potential 13| is modified by multiplying a factor 0.986 (1.365) for obtaining the 
desired value An ~ —0.3 (—9.0) MeV. Our result is very similar to that shown in Fig. 11 of 



Ref. [29J ; it reproduces the feature of stronger dependence of moment of inertia on pairing 
correlation for weakly bound systems. However, the free contribution j7o is rather small. 
In fact, the result J'^ > for the weak binding situation (An ~ —0.3 MeV) might seem 
strange at first sight because j7o in Eq. fl2B]) is always positive. Note that the Kruppa- 
BCS equation is solved for the calculation of instead of the ordinary BCS equation 
for J . Then the chemical potential obtained by the Kruppa-BCS equation is nearer 
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to the particle threshold, and the moment of inertia J becomes larger. The difference, 
j7(A(Kruppa-BCS))— J'(A(ordinary-BCS)), is larger than j7o, leading to the final result, 
> J . This clearly shows that even if JTq > its contribution is minor. In fact, for 
the stable (deeply bound) situation (An ~ —9.0 MeV) the difference between the Kruppa 
and ordinary calculations is negligible as is seen in Fig. |H1 We have done systematic 
calculations for many heavier nuclei 20| and this is generally the case. Namely, the 
neutron gas problem is not so harmful as far as the calculation of moment of inertia is 
concerned; subtraction of the free contribution in Eq. (126|) is not necessarily required. This 
is mainly because the free matrix elements of the operator in Eq. ( !27l) are small, which 
reflects that the momentum operators contained in the angular momentum {I = r x p) 
do not favor spatially-extended free wave functions in contrast to the spatial observable 
like the radius (r"^). 



III. SUMMARY 



We have extended the original Kruppa prescription for the single-particle level density 
to the calculation of one-body spatial observables like nuclear radius in the presence 
of pairing correlations. By simply subtracting the contribution of the free spectra, the 
effects of continuum states can be properly taken into account. The convergence property 
as increasing the basis size is carefully examined both numerically and analytically, and 
it is confirmed that reasonable convergence is attained for radius with manageable basis 
sizes except for the case of the halo-like situation with weak pairing correlation. 

The results of a few applications of the Kruppa method are presented for the basic 
quantities like nuclear radii, quadrupole moments, the neutron skin thickness, and the 
moment of inertia in very neutron rich nuclei. As for the moment of inertia, which is 
essentially the expectation value of the angular momentum operator at finite rotational 
frequency and is considered to be an example of operators consisting of both the coordinate 
and momentum, the effect of subtracting the free contributions is small. It has been 
found that the pairing anti-halo like effect naturally emerges, although the underlying 
mechanism is quite different. 

In this way, it has been shown that the so-called neutron gas problem is circumvented. 
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and the new Kruppa-BCS method gives rehable results similar to those of more sophisti- 
cated HFB calculations. Combined with the improved microscopic-macroscopic method 
developed in Ref . [lO| , we hope that the Kruppa prescription provides a simple, yet useful, 
new method for investigating nuclei far from stability. 

For the continuum single-particle level density, the Kruppa prescription in Eq. ([T]) has 
a sound basis . However, the replacement in Eq. (jS]) proposed in this work is 

missing such a solid foundation. It is an interesting future subject to clarify what is the 
meaning of this replacement, or on what kind of approximations it is based, from more 
fundamental many-body theory. 
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